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We study stability and structure of quark matters as a function of density in a frame- 
work of molecular dynamics (MD). Using appropriate effective interactions and the fric- 
tional cooling method, we search for the minimum energy of the system. Transition from 
confined to deconfined phase is observed at densities of 2 - 3p , where p is the nuclear 
matter saturation density. The uds matter becomes more stable than the charge-neutral 
Cj ' ud matter at 3p , but the udd matter is the most stable even at high density. 
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r'-J . 1. Introduction 
X 

More than twenty years ago Bodmer and later Witten put a hypothesis that the abso- 
lute ground state of strongly interacting matter is the deconfined state of quark matter 
consisting of an equal proportion of up, down, and strange quarks called strange quark 
matter even at zero temperature and zero pressure jl!2j . 

Recently some exotic baryon systems including strange quarks were discovered. For 
example, the penta-quark, called <9 + , consisting of 5 quarks (uudds) with an extraordinary 
small width ~25MeV, was discovered at SPring8 [3]. At KEK, a strongly bounded kaonic 
system, ppnK~ 1415 j . which was predicted from an AMD calculation [6 to form a high 
density system of 10po ; was discovered. From an astrophysical point of view, there is a 
possibility that the strange quark matter is realized at the center of neutron stars and/or 
quark stars [7 3 . At high baryon density such as in the core of these stars, quark gluon 
plasma (QGP) may exist and confined baryon and deconfined baryon states may coexist 
|8|9j . Searching for the energy region in which the hadron to quark transition occurs and 
the strange quark matter becomes stable is one of the most interesting topics in high 
energy hadron physics. 

In this paper we study the stability of quark matter at finite baryon density at zero 
temperature using a molecular dynamics method (MD). In this framework it is not nec- 
essary to assume whether the state of quarks is confined or deconfined and their mixed 
state can also be treated [TUj . 
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Figure 1. Snapshots for the udd matter after cooling. The white circles denote confined 
quarks and colored circles are deconfined quarks. 



2. Molecular dynamics for quark system 



We define the total wave function of a system \l/ with baryon number A as a direct 
product of single particle wave functions with a Gaussian wave packet in coordinate space 
and a state vector \i with a fixed spin orientation, flavor and color, 
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The antisymmetrization of the wave function is neglected, but treated phenomenologically 
by introducing a Pauli potential which acts as a repulsive force between quarks of the same 
color, flavor and spin orientation. The Hamiltonian is expressed as H = Ho + Vp an u— T S1 



- spur ; 



where T spur is a spurious zero point kinetic energy of clusters as described in [TT]. The 
original Hamiltonian H and Vp au ii are explicitly expressed as, 
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where f = |rj — Tj\ and A" = (Xi\^ a \Xi)- We use a s =1.25 (the QCD fine structure 
constant), iC=900MeV/fm (string tension), constituent quark masses m Uj( i=300MeV and 
m s =500MeV throughout this simulation which are typical values of parameters in quark 
models. The product of Gell-Mann matrices for different colors and for same colors 
are (E« = i ^i^j) = — I an d f j respectively However, since the antisymmetrization is ne- 
glected in our model, we multiply those values by 4 to effectively include the exchange term 
of color interaction. To introduce the nuclear force between white (colorless) baryons, we 
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Figure 2. Energy per quark 
for three kinds of matter. 
Dashed lines denote the result 
of "hadron-cooling" in which 
3 quarks always form a color- 
less hadron, while solid lines 
are those of ordinary "quark- 
cooling" without any constraint. 



employ the meson exchange potential Vmeson acting between quarks. The a and u meson- 
quark coupling constants are g aq =3A2 and 5^=9.0, the meson masses are /io-=300MeV 
and // li ,=782MeV, respectively. The parameters in the Pauli potential Vp au ii are deter- 
mined to reproduce the kinetic energy of the Fermi gas ^2j. In the present study we 
adopt C p =131MeV, g =l-8fm, p =120MeV for u and d quarks and C p =79MeV for s 
quarks. The size of the quark wave packet, L, is taken to be 0.43fm, whereas L eff =0.7fm 
is used in (y meS onj to take into account of the baryon size [TUj . 

3. Quark matter properties 

Here we investigate the ground state (zero temperature) of matter consisting of quarks. 
As an initial condition we distribute "hadrons" (made from 3 quarks with different colors) 
randomly in a cell with a periodic boundary condition. In the case of baryon density, e.g., 
of lpo; 90 quarks are put in a cell of a size 5.6 fm. To search for the minimum energy 
configuration of the system, we solve the cooling equations of motion with friction terms 
as in j!2j . The simulations are performed for u ud matter" (with same numbers of u and d 
quarks), "udd matter" (number of d quarks is twice that of u quarks) and ll uds matter" 
(with same numbers of u, d, and s quarks). 

Snapshots of the udd matter for different baryon densities are displayed in Fig. 1. All 
quarks are confined as hadrons at baryon density of lp . At 2.5p , a partly deconfined 
state can be seen, whereas almost all quarks are deconfined at 5p . 

The ground state energies of matter are shown with solid lines in Fig. 2. The dashed 
lines show results with a constraint that Rj of 3 quarks in each initial "hadron" should 
always be equal (to keep hadrons by hand). Here we call this calculation "hadron-cooling" , 
while the normal one we call "quark-cooling" . At lower density, the energy of the system 
by hadron-cooling agrees that by quark-cooling. This means that quarks are confined in 
hadrons as seen in Fig. 1. Note that the saturation of hadronic matter is realized around 
lpo for the first time in this kind of quark MD simulation. The energies of the ud matter 
and the udd matter are almost the same. This is because the Pauli potential is not fully 
effective for "hadrons" and p meson exchange is missing. At densities of 2 - 3po, transition 
to deconfined phase occurs. In our simulation, the decrease of the kinetic energy and the 
increase of the potential energy act oppositely during the deconfinement. The softening 
of the latter effect may be the origin of the transition at higher densities. 

Figure 3 shows comparison of energy among the ud matter with relativistic electrons, 
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Figure 3. Baryon density dependence 
of the energy per quark for three kinds 
of matter under charge-neutral condi- 
tion. 



the udd matter and the uds matter under charge neutral condition. Although the uds 
matter is the most unstable at low density, it becomes more stable than the ud matter as 
deconfmement occurs at 3po- The udd matter, however, is the most stable over all ranges 
in our present result. 



4. Summary 

Quark many-body systems were studied in a framework of MD. It has been found 
that hadron matter undergoes a phase transition to the quark matter as baryon density 
increases. At low density, the saturation for the ud matter is realized by introducing 
Vmeson- The udd and uds matter, however, also show the similar saturation because 
exactly the same meson exchange potentials are used among all quarks. To improve this, 
a flavor-dependent potential like p meson exchange potential is needed. The uds matter 
becomes more stable than the charge-neutral ud matter at 3po- But the udd matter is the 
most stable even at high density. We expect more reasonable results if p meson exchange 
potential is included in the future. 
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